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ABSTRACT 

We construct models for gamma-ray bursts where the emission comes from internal 
shocks in a relativistic wind with a highly non uniform distribution of the Lorentz 
factor. We follow the evolution of the wind using a very simplified approach where a 
large number of layers interact by direct collisions but where all pressure waves have 
been suppressed. We suppose that the magnetic field and the electron Lorentz factor 
reach large equipartition values in the shocks. Synchrotron photons emitted by the 
relativistic electrons have a typical energy in the gamma-ray range in the observer 
frame. Synthetic bursts are constructed as the sum of the contributions from all the 
internal elementary shocks and their temporal and spectral properties are compared 
to the observations. We reproduce the diversity of burst profiles, the "FRED" shape 
of individual pulses and the short time scale variability. Synthetic bursts also satisfy 
the duration-hardness relation and individual pulses are found to be narrower at high 
energy, in agreement with the observations. These results suggest that internal shocks 
in a relativistic wind may indeed be at the origin of gamma-ray bursts. A potential 
problem is however the relatively low efficiency of the dissipation process. If the rela- 
tivistic wind is powered by accretion from a disc to a stellar mass black hole it implies 
that a substantial fraction of the available energy is injected into the wind. 

Key words: gamma-rays: bursts; radiation mechanisms: non thermal; shock waves; 
accretion discs. 



1 INTRODUCTION 

Since 1991 the BATSE experiment on board the Comp- 
ton GRO satellite has observed more than 1800 gamma-ray 
bursts (hereafter GRBs). The burst distribution is isotropic 
over the sky but non homogeneous in distance (Fishman 
and Meegan, 1995 and references therein) which has been 
regarded as a strong indication that GRBs lie at cosmo- 
logical distances (Paczyhski 1991). However the possibility 
that GRBs belong to a large galactic halo (Hartmann et 
al. 1994) could not be excluded a priori. This long stand- 
ing controversy (see Nemiroff et al. 1995) about the burst 
distance scale may finally be solved by the recent observa- 
tions of transient optical counterparts for two GRBs. In the 
case of GRB 970228 (van Paradijs et al. 1997; Sahu et al. 
1997a) the point-like counterpart appears to be associated 
to an extended source which is probably a distant galaxy. 
The case of GRB 970508 is even more spectacular since the 
spectrum of the counterpart shows Fe II and Mg II lines 
due to absorbers on the line of sight at a redshift z — 0.835 
(Metzger et al. 1997). If the association of the visible source 
to the burst is confirmed, GRB 970508 must be a very dis- 
tant object at z > 0.835. If GRBs are placed at cosmologi- 
cal distances the LogiV - LogP (peak flux) curve and the 
value of < F/ Fmax > can be naturally interpreted in terms 



of cosmological effects and indicate that GRBs have typical 
redshifts in the range 0.3 - 1 (Piran 1992; Mao & Paczyhski 
1992; Fenimore et al. 1993). 

Modelling GRBs is a difficult task due to the extreme di- 
versity of burst profiles, the non thermal spectra and the lack 
of clear signature for the emission processes involved. Most 
cosmological models however share some common characte- 
ristics. The source, which must be able to release (between 
10 keV and 10 MeV) an energy _E 7 <; iiL- erg.sr -1 on a time 
scale of seconds, is generally supposed to be a stellar mass 
black hole accreting material from a disc. Such a configu- 
ration can result from the coalescence of two neutron stars 
(Eichler et al. 1989; Paczyhski 1991; Narayan, Paczyhski & 
Piran 1992) , the disruption of the neutron star in a neutron 
star - black hole binary (Narayan Paczyhski & Piran 1992; 
Mochkovitch et al. 1993) or the collapse of a massive star 
(Woosley 1993). The power emitted by cosmological GRBs 
is orders of magnitude larger than the Eddington luminosity 
and cannot come directly from the disc surface. The released 
energy instead drives a wind which has to become relativis- 
tic both to produce gamma-rays and to avoid photon-photon 
annihilation along the line of sight (Baring 1995; Sari & 
Piran 1997a). The Lorentz factor P must reach values of 
10 2 - 10 3 which limits the allowed amount of baryonic pol- 



2 F. Daigne and R.Mochkovitch 



lution in the flow to a very low level. A few mechanisms 
which could possibly achieve such a severe constraint have 
been proposed: (i) magnetically driven outflow originating 
from the disc or powered by the Blandford-Znajek (1977) 
effect (Thompson 1994; Meszaros & Rees 1997a), ii) recon- 
nection of magnetic field lines in the disc corona (Narayan, 
Paczynski & Piran 1992) or (iii) neutrino-antineutrino an- 
nihilation in a funnel along the rotation axis of the system 
(Meszaros & Rees 1992; Mochkovitch et al. 1993, 1995). It 
is supposed in (i) and (ii) that the field has reached huge 
values B £ 10 15 G. Concerning (iii) Ruffert et al. (1997) 
have shown recently that vv annihilation does not provide 
enough energy to account for cosmological GRBs except may 
be for very massive discs, as those which could result from 
the collapse of a massive star. 

The energy initially stored in kinetic form within the 
relativistic wind must then be converted into gamma-rays. 
This can be done during the deceleration of the wind resul- 
ting from its interaction with the interstellar medium (Rees 
& Meszaros 1992; Meszaros & Rees 1993) or a dense radia- 
tion field (Shemi 1994). In the first case the emission comes 
from electrons accelerated in the forward and reverse shocks 
which then radiate synchrotron and inverse Compton pho- 
tons in a magnetic field frozen in the wind or which has 
come to equipartition with the shocked material (Meszaros, 
Laguna & Rees 1993). In the second case, photons of the 
ambient radiation field (in the central regions of globular 
clusters or active galactic nuclei) interact with the electrons 
(Shemi 1994; Shaviv & Dar 1995) or ions (Shaviv & Dar 
1996) of the wind and undergo a boost in energy by a fac- 
tor r 2 . With F 10 3 , optical photons can be shifted to 
the gamma-ray range. In these two models the duration of 
the burst is t\, ~ ^pf- where r^cc is the deceleration ra- 
dius of the wind. Another possibility suggested by Rees & 
Meszaros (1994) consists to suppose that the Lorentz fac- 
tor in the wind is variable so that successive shells can have 
large relative velocities leading to the formation of internal 
shocks. The energy which is dissipated in these shocks can 
then be radiated as gamma-rays via the production of pi- 
ons in proton-proton collisions (Paczynski & Xu 1994) or by 
synchrotron (or inverse Compton) emission of accelerated 
electrons. There are several important differences between 
the deceleration and internal shock models. In deceleration 
models, the duration and time profile of the burst depend 
on r, on the value of the deceleration radius and on the 
structure of the emitting shell. In internal shock models, 
the duration of the burst is directly related to the duration 
of energy injection at the source and the time profile is es- 
sentially determined by the variations of the Lorentz factor 
(Sari & Piran 1997a). 

Burst profiles and spectra have already been obtained 
(sometimes in a rather detailed manner) in the case of de- 
celeration models (Fenimore, Madras & Nayakshin 1997; 
Panaitescu et al. 1997; Panaitescu & Meszaros 1997) and the 
purpose of this work is to present the same kind of quantita- 
tive analysis for internal shock models. In Sect. 2 we describe 
our method to follow the evolution of the relativistic wind 
and we discuss the emission processes; Sect. 3 and 4 respec- 
tively deal with the temporal and spectral properties of our 
synthetic burst models and Sect. 5 is the conclusion. 



2 A SIMPLE MODEL OF THE RELATIVISTIC 
WIND 

2.1 Description of the model 

We do not discuss in this paper the nature of the source 
(coalescence of two neutron stars, neutron star - black hole 
binary, collapse of a massive star or something else) which is 
initially responsible for the energy release. We also suppose 
that a relativistic wind carrying the energy has emerged 
from the source, with an average Lorentz factor f of a few 
hundreds. 

To study the evolution of the relativistic wind we have 
developed a very simple model where a succession of layers 
are emitted every 2 ms with a varying Lorentz factor, during 
a total time t w . The mass of the layers is proportional to 1/r 
so that the energy injection rate is constant. We follow the 
layers as the wind expands and when a rapid layer (of mass 
mi and Lorentz factor Ti) catches up with a slower one 
(ni2, T2 < Ti) they collide and merge to form a single shell 
of resulting Lorentz factor r r . If the dissipated energy 



e = m\<?Y\ + m,2C 2 F2 



(mi + ?Tl2)c r r 



(1) 



can be radiated in the gamma-ray range on a time scale 
shorter than the shell expansion time (in the comoving 
frame) 



r 



(2) 



where r is the shell radius, the burst profile will be made by 
a succession of elementary contributions of duration (in the 
observer frame) 



At: 



2cr r 2 



(3) 



We estimate T r by considering that most of the energy avail- 
able in the collision has been already released when the less 
massive of the two layers has swept up a mass comparable to 
its own mass in the other layer (internal shocks being only 
mildly relativistic). Then 



and 

e = mc (Vi + T 2 - 2r r 



m = min(mi, 7712) 



(4) 



(5) 



After the complete redistribution of momentum and energy 
the Lorentz factor of the merged layer finally becomes 



miTi + m,2r 2 



riP 



miF2 + m^Ti 



(6) 



Two more conditions have to be satisfied for an elementary 
shock to be produced and observed: i) the two layers must 
collide at a relative velocity larger than the local sound speed 
and ii) the wind must be transparent to the emitted photons. 
The relative velocity of the two layers is given by 



c 



r 2 + r 2 



(7) 



with Ti > P2 and IY2 3> 1. We have adopted a sound speed 
v s /c — 0.1 but we have checked that other choices make little 
differences in the results since the main contribution to the 
burst comes from shocks corresponding to large differences 

of the Lorentz factor <J 2. 
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The transparency of the wind to the emitted photons 
has been computed in the following way: in each collision 
between layers a "photon shell" is generated which then 
catches up with all the layers ahead of it. It travels through 
a total optical depth 



rrii 
47rr 2 



(8) 



where kt is the Thomson opacity and n the radius of layer 
i when it is reached by the photon shell. The sum is over 
all indices i larger than i m , which corresponds to the two 
merged layers. For the average Lorentz factor f J> 100 used 
here it appears that the wind is transparent to the emitted 
photons except for a few early collisions. 

The arrival time of each of the elementary contributions 
from internal shocks is calculated relatively to a signal which 
would have travelled at the speed of light from the source 
to the observer. It is given by 



ta — te 



(9) 



where t e is the emission time and r the distance to the source 
of the two merged layers. The evolution of the system is 
followed until all the layers are ordered with F decreasing 
from the front to the back of the wind. The efficiency of the 
dissipation process can then be obtained as 



/d = 



(10) 



where the e s are the energies released in each of the internal 
elementary shocks and the m,, Ti are the initial masses and 
Lorentz factors of the layers. 

This very simple approach is naturally very crude be- 
cause it neglects all pressure waves propagating throughout 
the layers. Nevertheless, we expect that it can still capture 
the basic features of the real process. As a first example 
we have represented in Fig.l the evolution of the Lorentz 
factor at different times when the initial distribution of F 
consists of 5000 layers (i.e. t w = 10 s) with T(n) = 400 for 
n = 1 to 4000 and T(n) = 100 for n = 4001 to 5000 (n = 1 
for the last emitted layer). Such a wind is then divided be- 
tween a "slow" part which is emitted first and a rapid part 
which will progressively collide against the slow part. The 
total masses injected in the slow and rapid parts are equal. 
After the first collision has occurred a separation layer of 
increasing mass is formed between the slow and rapid parts. 
We have plotted in Fig. 2 the values of t e , F r , At and e for 
the elementary shocks as a function of arrival time t a . The 
evolution of the shell system is essentially completed after a 
time t e ~ ivarf 2 , where t var is the characteristic time scale 
for the variations of the Lorentz factor. As expected, the 
plots of T r , At and e have two branches, corresponding to 
collisions taking place on both sides of the separation layer. 
They mimic the forward and reverse shocks which, in our 
simple model, are located on a single surface. 

The efficiency /d of the dissipation process is rather 
low, 10% in this specific case, typically less than 15%. This 
is a severe problem, which can become even worse if only a 
fraction of the dissipated energy is radiated in the gamma- 
ray range. This issue is examined in the next section. 



Figure 1. Distribution of the Lorentz factor in the wind at differ- 
ent times in lagrangian coordinates m/rritot (m-tot being the total 
mass of the wind). The initial distribution (t = 0) consists of 5000 
layers with T = 100 in the 1000 first emitted layers and T = 400 
in the rest of the wind. Since the mass of the layers is propor- 
tional to 1/r the slow and rapid parts contain an equal mass. At 
t = 3.4 10 4 s a forward and a reverse shock are propagating into 
the wind. In our simple model they are located on a single eule- 
rian shell of increasing mass. At t = 1.7 10 5 s the forward shock 
has already crossed the whole slow part and at t = 5.4 10 5 s the 
reverse shock reaches the end of the rapid part. The dashed line 
shows the initial distribution of the Lorentz factor. 

Figure 2. Physical parameters of the elementary shocks. Upper 
panel: (thick line) emission time t c and (thin line) duration At 
of the signal (in the observer frame) both as a function of arrival 
time t a . Lower panel: (thick line) energy e dissipated in the shock 
in arbitrary unit and (thin line) Lorentz factor T r at the shock 
location. 



2.2 Emission processes 

The process by which the dissipated energy is finally ra- 
diated depends on the energy distribution of protons and 
electrons in the shocked material and on the values of the 
comoving density and magnetic field. The average energy 
which is dissipated per proton in a shock between two lay- 
ers of equal mass is given by 



e = (r int - l)r 



with Tin 



® 



1/2 



1/2' 



(11) 



where Tint is the Lorentz factor for internal motions in the 
shocked material. With Fi/F2 = 4, which corresponds to a 
mildly relativistic shock (v le i/c = 0.88), e > 200 MeV. 

A mechanism which could directly extract the energy 
from the protons has been proposed by Paczynski & Xu 
(1994) since the value of e is large enough for pions produc- 
tion in pp collisions. The pions then decay with the emission 
of gamma-rays. The global efficiency of this mechanism is 
low, of the order of 10~ 3 for V 1 /T 2 = 4. 

Another possibility, considered by Rees & Meszaros 
(1994), Papathanassiou & Meszaros (1996) and Sari & Piran 
(1997a), consists to suppose that the electrons have come 
into (at least partial) equipartition with the protons. If a 
fraction a c of the dissipated energy goes to the electrons 
their characteristic Lorentz factor will be 

(12) 



m c c^ 



which, for ct e = 1/3 (corresponding to a complete equiparti- 
tion between protons, electrons and the magnetic field) and 
ri/r2 = 4 yields F c ~ 150. The equipartition magnetic field 
is given by 



Boq ^ (87raBne) 1/2 



(13) 



where qb Si 1 and n is the comoving proton number density 



M 



E 



47rr 2 rm p c 



(14) 



Assuming E — 10 52 erg.s 1 , V = 300, t var = 1 s and qb = 
1/3 the equipartition magnetic field at a radius r ~ ct var f 2 
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where most of the collisions take place is B cq ~ (10 2 — 10 3 ) G 
depending on the ratio ri/IY 

Synchrotron emission by the accelerated electrons in the 
magnetic field occurs at a typical energy (in the observer 
frame) 

which corresponds to the UV range for T c = 100. Gamma- 
rays can be produced by inverse Compton scattering on the 
synchrotron photons. Then 



E IC ~ E Byn T 2 e = 500 



1300/ U000G/ 1 100 J 



keV, (16) 



and the fraction of the total power which is radiated by the 
inverse Compton process is 

aic = , (17) 



i + r*r c 2 



where 



47rr 2 



(18) 



is the optical depth of the shell of mass M* and radius r* 
which contains the relativistic electrons. An estimate of M* 
is 



Mi, = 



1 + T*H 



M s hock , 



where 



*^ = 6 (l^) 



B 



\ 1000G 



(19) 



(20) 



is the synchrotron time of the relativistic electrons and 
A^shock the mass flow rate across the shock, both in the co- 
moving frame of the shocked material. Since the shock moves 
with a Lorentz factor r s ~ f , M shoc k can be approximated 
by 



Mshock 



M 

7 



(21) 



From (18), (19) and (21) we obtain an implicit expression 
for r*r 2 



47rr 2 (l + r*r 2 ) 



(22) 



With r* ~ ct var r 2 , relations (20) and (21) for t syn and 
Mhock, (13) and (14) for B cq and n, Eq. 22 gives (for 

KT 



: 0.2 cm 2 .g~ lN 



r,r 2 (l + r,r 2 )~0.3^ , (23) 

QB 

which, for a c ~ qb yields r*r 2 ~ 0.24 and aic — 0.19. A 
larger fraction of the dissipated energy can be converted to 
gamma-rays if the magnetic field does not reach equiparti- 
tion. With for example ^ = 100, t*T 2 rises to 5 and aic 
to 0.83. Smaller values of qb could still increase aic but the 
energy of the inverse Compton photons would then become 
too small. 

Conversely, if the electron Lorentz factor is large enough 
gamma-rays can be directly produced by synchrotron emis- 
sion. This will be the case if, instead of (12) one uses the 
expression of T c given by Bykov & Meszaros (1996) who 



consider the scattering of electrons by turbulent magnetic 
field fluctuations. They get 

l/(3-/i) 



T v 



(24) 



where a M is the fraction of the dissipated energy which goes 
into magnetic fluctuations, £ the fraction of the electrons 
which are accelerated and \i the index of the fluctuation 
spectrum. With 1.5 < [i < 2, a M = 0.1 - 1, C ~ 10~ 3 and 
e/m e c 2 ~ 500 (for Ti/r2 = 4) values of T c in the range 
10 3 - 10 4 can be obtained. 

A fraction of the synchrotron photons will be shifted 
to even higher energy by inverse Compton scattering which 
now occurs in the limit where 



r p° 

1 c-C^svr 



~ 33 



B 



1000G 



(l0 4 ) 



(25) 



is large, -E° yn = E syn /r T being the synchrotron energy in 
the comoving frame. The inverse Compton photons have an 
energy 

2 = 5 



^ic 



r o m c 



GeV 



VIO 4 / 

in the comoving frame and carry a fraction 



aic 



(26) 



(27) 



1 + r.Yl/w ' 

of the dissipated energy. The optical depth r* has to be 
computed with the Klein-Nishina cross-section which, in the 
limit w 3> 1, gives 



KT-Mshockis 



47rr 2 (l + r*r 2 /™) 
and therefore 



(I) t^ 1 ^)], 



(28) 



r*r c 2 



i + 



8 10 _4 [1 + ln(2w)] 



E 



10 52 erg.s -1 



\T7J 



300 J llOOOcJ llO 4 ) 



(29) 



The very large exponents which appears in (29) show that a 
small variation of the parameters can induce large changes 
in the relative importance of the synchrotron and inverse 
Compton processes. In practice, we find that aic is general- 
ly small in the early part of a burst but increases at later 
times essentially due to the reduction of the equipartition 
magnetic field in shocks at large distances from the source 
(see Sect. 3.1 below). 

To be efficient the emission process must also occur on 
a time scale i em shorter than the shell expansion time t ex 
(Eq. 2). Using r ~ ct var r 2 , this condition becomes 

U 



< rt v 



(30) 



1 + Qic 

where Qic = r*r 2 (resp. tvX 2 /™) for w <§C 1 (resp. w >• 1). 
With expressions (11), (13), (14) and (20) above, (30) can 
be written 



2 10- 4 a B 1 (r int - I)" 1 
't \ I f 



10 52 erg.s -1 



fe) (m &r 



(31) 
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It is more easily satisfied when gamma-rays directly come 
from synchrotron emission since then r e ~ 10 4 . However, 
even for T c ~ 10 2 the emission time remains smaller than 
the expansion time as long as the shocks are sufficiently 
strong (T1/T2 k, 2). 

Some results of a model where GRBs are produced by 
the inverse Compton process have been already presented 
elsewhere (Mochkovitch & Fuchs, 1996). In this study we 
limit ourselves to synchrotron emission models except when 
we discuss the optical properties which strongly differ be- 
tween the two cases. In synchrotron emission models the 
total efficiency for the conversion of wind kinetic energy to 
gamma-rays below a few MeV is given by 



/tot = fd x a e (l — aic) . 



(32) 



where a c is the fraction of the dissipated energy which is 
transferred to the electrons. According to Bykov & Meszaros 
(1996) Q c is comparable to the fraction bm ~ 0.1 — 1 of the 
energy which is initially injected into magnetic fluctuations. 
The total efficiency therefore does not exceed a few per- 
cent which imposes severe constraints on the energy source 
or/and system geometry. 

The energy available from disc accretion to a black hole 



E ~ X -Muc 2 = 3 10 53 ( ^ 
6 \M e 



erg 



(33) 



where Md is the disc mass, can be less than 10 53 erg for 
the coalescence of two neutron stars since numerical simu- 
lations indicate that Mo ~ 0.2 — 0.3 Mq (Rasio & Shapiro 
1994; Davies et al. 1994; Ruffert, Janka & Schaefer 1996). 
The conversion of disc gravitational energy into wind kinetic 
energy should then be very efficient in order to account for 
cosmological GRBs. Naturally, if the wind is beamed in a 
solid angle 80, along the system axis the energy requirement 
is smaller by a factor |^ but one has now to face a statistical 

J 27T 

problem since even optimistic estimates of the neutron star 
merging rate are not considerably larger than the burst rate 
(Phinney 1991; Narayan, Piran & Shemi 1991; Tutukov & 
Yungelson 1993; Lipounov et al. 1995). The situation is less 
critical if the disc results from the disruption of a neutron 
star by a black hole or the collapse of a massive star because 
the mass of the disc can be Mo <; 1 M Q (and even Mo ii 10 
M ). 



3 TEMPORAL PROPERTIES 
3.1 Burst profiles 

We first study the temporal properties of our burst models 
when the Lorentz factor in the wind has the simple shape 
shown in Fig.l. We inject an energy E — i2— erg.s -1 .sr -1 
which, for an efficiency of conversion into gamma-rays of 
a few percent, yields a total energy E ~ erg.sr -1 for 
a burst lasting a few seconds. We obtain the burst profiles 
(in number of photons per second between 50 and 300 keV, 
which corresponds to BATSE bands 2 and 3) by adding the 
contributions of all the internal elementary shocks which 
occur during the expansion of the wind as explained in Sect. 



2. For each shock we compute a B 



1 - aic (Eq. 27) 



Figure 3. Physical parameters governing the emission mecha- 
nism: (a) fraction o S yn = 1 — ctic °f the energy which is radiated 
by the synchrotron process; (b) Lorentz factor T e of the relativis- 
tic electrons; (c) cquipartition magnetic field B cq ; (d) synchrotron 
energy E Byn . In (a), (b) and (d) the thick lines correspond to T e 
given by Eq. 24 and the thin lines to r o = 10 4 . 



Figure 4. Burst profiles for the initial distribution of the Lorentz 
factor shown in Fig.l. The count rate C2+3 (in arbitrary unit) 
is given in the interval 50 - 300 keV, corresponding to BATSE 
bands 2 and 3; (a) profile obtained with Eq. 24 for the electron 
Lorentz factor; (b) same as (a) with C2+3 in logarithmic scale 
which illustrates the exponential decay after maximum; (c) profile 
obtained with a constant r o = 10 4 ; (d) same as (c) with C2+3 in 
logarithmic scale. 

13 with cvb = 1/3), the electron Lorentz factor T e (Eq. 24 
with = 1000 and fi = 1.75) and the synchrotron energy 
E syn (Eq. 15). These quantities are represented in Fig. 3 as a 
function of arrival time £ a . As in Fig. 2 there are two branches 
corresponding to the forward and reverse shocks. Typical 
values are 10 2 - 10 4 G for the magnetic field, 2 - 20 x 10 3 
for the electron Lorentz factor and 10 keV - 1 MeV for the 
synchrotron energy. 

The photons produced in the elementary shocks are dis- 
tributed according to a synchrotron spectrum 



ME) rr g 
dE 



Es- 



E 



(34) 



where aic is the fraction of the energy which goes to inverse 
Compton photons, the equipartition magnetic field B cq (Eq. 



with x = 2/3 for E < E syn , 2 < x < 3 for E > E syn and 
where e is the energy which is dissipated in the shock (Eq. 
5). The resulting burst profile is shown in Fig. 4a for a high e- 
nergy index x = 2.5. Cosmological effects (time dilation and 
redshift) have been taken into account, assuming that the 
burst is located at z = 0.5. The duration tgo which, accor- 
ding to the definition used for the BATSE data is the time 
during which 90% of the total fluence is received (excluding 
the first and last 5%) is 10.04 s, very similar to the duration 
of wind emission t w . The decay after maximum is close to an 
exponential as can be seen in the plot of the logarithm of the 
count rate versus time in Fig. 4b where there is a quasi linear 
decline after maximum between t ~ 12 s and t ~ 20 s. How- 
ever the rise time in our profile is not much shorter than the 
decay time and the burst is therefore not very dissymmetric 
if we except the low intensity exponential tail after t = 15 
s. We found that a profile much closer to the characteris- 
tic FRED (Fast Rise Exponential Decay) shape observed in 
many bursts (or in individual pulses inside a complex burst) 
can be obtained by assuming that the electron Lorentz factor 
varies more slowly with e (the dissipated energy per proton) 
than a power- law of index 1/(3 — n) (Eq. 24). This would be 
the case if the fraction ( of accelerated electrons increased 
with e. Such a behavior is observed in simulations of col- 
lisionless non relativistic shocks (Bykov, private communi- 
cation) and we have supposed that it remains valid in the 
relativistic limit. Moreover, we have made the simple choice 
C oc e which leads to a characteristic Lorentz factor for the 
electrons which is independent of e. The synchrotron energy 
in the elementary shocks has been represented in Fig. 3d for 
a constant r o = 10 4 and the corresponding burst profile is 
shown in Fig. 4c. It has now a typical FRED shape with a 



6 F. Daigne and R.Mochkovitch 



Figure 5. Full line: ratio of the decay time to the rise time of 
burst profiles obtained with an initial distribution of the Lorentz 
factor homothetic to that shown in Fig.l but with different du- 
rations. For tt)o < 1.65 s the profiles decay faster than they rise. 
Dashed line: the discontinuity in the initial distribution of the 
Lorentz factor has been replaced by a smoother transition (Eq. 
36). Now, only bursts with tgo < 0.65 s decay faster than they 



Figure 6. Burst profiles for three initial distributions of the 
Lorentz factor in the wind. In all three cases a rapid compo- 
nent with r = 400 is decelerated by a series of slower layers. The 
masses in the rapid component and slower layers are compara- 
ble, (a) Relatively simple profile with four layers which produce 
four intensity pulses two of which partially overlap; (b) and (c) 
more complex profiles with 15 slow layers. Notice that the distri- 
butions of the Lorentz factor in (b) and (c) are homothetic, (c) 
being simply ten times longer than (b). 

Figure 7. Same as Fig. 6 with a random fluctuation (of maximum 
amplitude ± 20%) added to the average value T = 400 of the 
rapid component of the wind. The resulting profiles now exhibit 
variability on a short time scale. 

ratio of the decay time to the rise time — = 3.4 where T r 



and Td are defined respectively by 



tr, 



t 5% 



(35) 



where i max is the time of maximum count rate and t 5 % (resp. 
^95%) the time when 5% (resp. 95%) of the total fluence 
has been received. We have tested with our model the ten- 
dency for short bursts (or short pulses within a complex 
burst) to become more symmetric (Norris et al. 1996). We 
have computed the profiles obtained when the duration of 
wind emission t w is varied while the initial distribution of 
the Lorentz factor remains homothetic to a given shape for 
which we choose the one represented in Fig.l where the wind 
consists of a slow part (F = 100) followed by a rapid part 
(r = 400) both containing the same mass. The results (for 
z — 0.5) are shown in Fig.5 where we have plotted the ratio 
as a function of tgo- It appears that ^f- decreases from 
about 3 when tg ~ 10 s to about 0.3 when tg <, 0.5 s 
with — = 1 for tgo ~ 1.65 s. We therefore reproduce the 
observed behavior but the effect is even exagerated for the 
shortest pulses which decay faster than they rise. We believe 
that this might be a consequence of the crude assumptions 
made in our simple model and we expect that more detailed 
hydrodynamical simulations (Daigne and Mochkovitch, in 
preparation) could help to improve the profiles. Already, we 
found that better results can be obtained (with ^ = 1 for 
tgo — 0.65 s) if the discontinuity between the two extreme 
values of F is replaced by a smoother transition of the form 



r(t/t w ) = 250 + 150 cos[2.57r(i/i w - 0.6)] , 



(36) 



for 0.6 < t/t v < 1 and T = 400 for t/t w < 0.6. However, 
the efficiency /d of the dissipation process is then reduced 
by nearly a factor of 2. 

Norris et al. (1996) have shown that in most cases com- 
plex bursts can be analyzed in terms of a series of (possi- 
bly overlapping) pulses. In the same way we build complex 
bursts with our model by the addition of intensity pulses 
formed in the deceleration of rapid parts of the wind by 
slower ones which were emitted previously. We suppose that 



the initial distribution of the Lorentz factor is made of a 
rapid component (with an average value of Y of a few hun- 
dreds) and of some slower layers (with Y ~ 100). The total 
mass in the slow layers has to be comparable to the mass 
in the rapid component in order to keep the efficiency at a 
reasonable level. A few examples of synthetic profiles are pre- 
sented in Fig. 6. It can be seen that a great diversity of burst 
shapes can be obtained if the distribution of the Lorentz 
factor in the wind varies from one event to the other. 

3.2 Short timescale variability 

The profiles shown in Fig. 4 and 6 have a satisfactory gene- 
ral shape but do not exhibit any variability on a short time 
scale. Rapid temporal variations cannot result from small 
irregularities of the emitting surface since photons coming 
from many different regions in space and time are received 
at a same time by the observer which leads to a loss of cohe- 
rence of the temporal variations (Woods & Loeb 1995; Sari 
& Piran 1997b, c). Instead, one has to consider again that the 
Lorentz factor itself can fluctuate at the millisecond level. 
This is not unrealistic if the flow at the origin of the wind 
is very irregular and turbulent since one millisecond corres- 
ponds to the typical dynamical time scale of a disc orbiting 
a stellar mass black hole. Nothing being known about the 
temporal spectrum of the fluctuations we have simply added 
a random fluctuation to the Lorentz factor of each of the lay- 
ers initially injected (every 2 ms) in the wind. The adopted 
amplitude for these fluctuations is 10 - 20% of the average 
value of T. The resulting profiles represented in Fig. 7 now 
show the rapid temporal variations which are seen in most 
observed bursts. 



4 SPECTRAL PROPERTIES 
4.1 Burst spectrum 

The overall burst spectrum is the sum of all the elemen- 
tary contributions (Eq. 34) from the internal shocks. The 
spectrum corresponding to the profile shown in Fig. 7a is 
represented in Fig. 8 (again a cosmological redshift z = 0.5 
has been assumed). Its shape can be easily understood: let 
and El^f be the minimum and maximum of the syn- 
chrotron energy for the whole set of elementary shocks. As 
long as E is smaller (resp. larger) than -E™' n n (resp. K^f) 
the number of photons n(E) per unit energy interval is a 
power-law of index —2/3 (resp. —2.5) and in the interme- 
diate region the current index evolves from —2/3 to —2.5 
as E goes beyond the value of E Byn for a growing number 
of elementary shocks. In the four BATSE bands, between 
20 keV and a few MeV, the spectrum can be well described 
with Band's formula (Band et al. 1993) 



i(E) = A ( 



n(E) = A 



E 



100 keV 

(a - 0)Eo 



exp 



100 keV 

for (a - P)E < E . 



( E ) 
exp(/3 - a) ( 



for (a - (3)Eo > E , 

r 



E 



100 keV, 



(37) 



The parameters a, (3 and Eq have been adjusted to ob- 
tain the best possible fit of the spectrum in Fig. 8. We get 
a = —1.33, P — —2.31 and -Bo = 544 keV, in agreement 
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Figure 8. Spectrum of the burst corresponding to the profile of 
Fig. 7a. The number of photons per energy interval n(E) and the 
product E 2 n(E) are shown in arbitrary units. The dashed line is 
a fit of the spectrum with Band's formula in the interval 10 keV 
— 10 MeV. The product E 2 n(E) is maximum at the peak energy 
E p = 365 keV. 

with typical values found in observed bursts. The product 
E 2 n(E) is also shown in Fig. 8. It is maximum at the peak 
energy E p = 365 keV where the bulk of the emission takes 
place. 



Figure 9. Duration-hardness ratio (HR32) relations for a simple 
(one pulse) burst and a complex burst with five pulses. The rela- 
tions are shown for three values of the high energy index x = 2, 
2.5 and 3 of the elementary synchrotron spectrum and two red- 
shifts z = 0.3 and 1. It can be seen that the effect of the redshift 
is negligible. 

Figure 10. Evolution of hardness with time: the instantaneous 
value of the peak energy E p (thick line) is represented together 
with the burst profile (thin line). Upper panel: for the simple one 
pulse burst of Fig.4c. Lower panel: for the more complex burst of 
Fig.6a. 



4.2 Duration-hardness ratio relation 

Shorter bursts are expected to be harder in our model. Inter- 
nal shocks are formed at an approximate radius r ~ ct var f 2 
and are therefore closer to the source if the burst evolves on 
a short time scale. Assuming that the injected power E is 
independent of the duration of wind injection t w the equipar- 
tition magnetic field is stronger and the synchrotron energy 
is larger in the dissipation region for shorter bursts. To ob- 
tain the duration-hardness relation we compute the hardness 
ratio HR32 (defined as the ratio of the number of photons re- 
ceived in BATSE band 3 to that in band 2) for bursts with 
homothetic initial distributions of the Lorentz factor but 
different values of t w - For a given duration tgo the hardness 
ratio is a function of the high energy index of the elemen- 
tary spectrum (Eq. 34) (the bursts becoming softer when x 
increases from 2 to 3), of the detailed history of the Lorentz 
factor and of the cosmological redshift z. Complex bursts 
tend to be harder because the distribution of the Lorentz 
factor varies on a shorter time scale at given t w and dissi- 
pation therefore begins earlier than in more regular bursts 
of same duration. Finally cosmological effects shift bursts in 
the HR32 - tgo diagram downward (redshift) and to the right 
(time dilation). However the shift is nearly colinear to the 
duration-hardness relation which therefore remains practi- 
cally unchanged at different z. This is illustrated in Fig. 9 
where the duration-hardness relation has been represented 
for several values of x and z and for a simple (single pulse) 
and a complex burst. In agreement with the observations 
(Kouveliotou et al. 1993; Dezalay et al. 1996) a transition 
occurs at tgo ~ 2 s. The shortest bursts reach a limit 

300 1 / 3 — 100 1//3 
HR32 ~ xooi/3 _ 50 i/3 ~ 2A > ( 38 ) 

when E p is large enough for both BATSE bands 2 and 3 to lie 
in the region where n(E) oc E^ 2 '' 3 . The longer bursts tend 
to various limiting values of the hardness ratio depending 
on the choice made for the high energy index x. 

4.3 Spectral evolution 

4-3.1 Instantaneous hardness 

The spectral evolution of GRBs shows a few trends which 
are followed by a majority of bursts but also suffer some ex- 
ceptions and are therefore not universal (Bhat et al. 1994; 
Ford et al. 1995). First, spectral hardness and count rate ap- 
pear to be correlated. Within intensity pulses both increase 
and decrease together, the hardness usually preceding the 
count rate. Another trend is a global hard-to-soft evolution 
over the course of the burst outside intensity pulses. Finally, 



later pulses tend to be softer than earlier pulses, even if they 
have a greater intensity. 

We have compared the spectral evolution of our syn- 
thetic burst models to these observational results. The hard- 
ness can be obtained as a function of time through the es- 
timation of the instantaneous value of E p , the energy of 
the peak of E 2 n(E). We first considered the simple burst 
of Fig. 4c which has a characteristic FRED shape and a du- 
ration tgo — 10.23 s. The hardness and count rate evolve 
similarly (see Fig. 10) but their maxima are separated by 
a time lag of 0.89 s, the maximum of E p occurring before 
that of the count rate. In complex bursts the hardness in- 
creases during intensity pulses and also precedes the count 
rate (Fig. 10). 

The correlation between spectral hardness and count 
rate is therefore correctly reproduced in our models. We en- 
countered more difficulties with the global hard-to-soft evo- 
lution which was observed in 70% of the sample of bright 
long bursts studied by Ford et al. (1995). In synthetic bursts 
(Fig. 10) it is present as long as the profiles remain relatively 
simple (i.e. dominated by one main pulse or made of just a 
few pulses) while in complex bursts with many pulses only 
the correlation between hardness and count rate is clearly 
visible. Also, the hardness of successive pulses remains es- 
sentially correlated to the intensity instead of decreasing like 
in 50% of the Ford et al. (1995) sample. 



4-3.2 Pulse shape as a function of energy 

When observed in spectral bands of increasing energy, pulses 
become narrower as shown by Norris et al. (1996) who ana- 
lyzed a large number of pulses in the four BATSE bands. 
They found that their (half maximum) width can be well 
represented by a power-law 



W(E) oc E~ 



(39) 



The same relation was obtained by Fenimore et al. (1995) 
who used the autocorrelation of averaged burst profiles in- 
stead of individual pulses. In our synthetic burst models the 
pulse width also decreases at high energy. Figure 11 repre- 
sents a single pulse burst in the four BATSE bands. The 
width can be fitted by power-laws such (39) but not with 
a unique exponent p = —0.4 for all the pulses. We indeed 
get p ~ —0.4 for pulses of 2 - 10 s but for pulses of 0.1 - 
1 s, p <; —0.2. Norris et al. (1996) obtained p = —0.4 as an 
average on a collection of pulses with duration ranging from 
0.1 to 10 s but do not provide the value of p for the shorter 
and longer pulses separately which does not allow a detailed 
comparison with our theoretical results. 
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Figure 11. Upper panel: normalized profiles for the burst of 
Fig. 4c now represented in all four BATSE bands. Lower panel: 
half maximum widths of the profiles as a function of energy. The 
fit of the model results (dashed line) has a slope p = —0.39. 



4.4 The E p - fluence relation 

Liang & Kargatis (1996) discovered in a sample of 37 BATSE 
bursts an exponential dependence of the peak energy E p on 
photon fluence F pb during the declining part of intensity 
pulses i.e. 

E p oc cxp(-aFph) . (40) 

The photon fluence is defined as the integral of the photon 
flux from the beginning of the burst and a is the slope of the 
LogiSp - Fph relation. In complex bursts the slope stays ap- 
proximately constant from one pulse to another. Synthetic 
bursts also follow a relation such (40) but the slope in suc- 
cessive pulses can somewhat vary, especially if they have 
very different intensity or duration (see Fig. 12). 

4.5 X-ray and optical counterparts 

The X-ray and optical counterparts recently discovered in 
two GRBs have been interpreted in the context of cosmo- 
logical models as the emission coming from a rclativistic 
shell expanding in the interstellar medium (Wijers, Rees & 
Meszaros 1997; Vietri 1997; Waxman, 1997a,b). If however 
GRBs are produced by internal shocks in a wind, X-ray to 
optical photons should also be emitted together with the 
gamma-rays before the afterglow resulting from the interac- 
tion with the interstellar medium. 

To compute these early counterparts we consider a ty- 
pical burst where an energy E — 2 *° erg.sr^ 1 has been 
injected into the wind with a 5% efficiency for the conver- 
sion to gamma-rays between 50 keV and 300 keV. For a 
GRB located at 2 Gpc (z ~ 0.5) the observed fluence in 
BATSE bands 2 and 3 will be F-, ~ 2 10~ 6 erg.cm" 2 . This 
value of Fj is then used to normalize a synthetic spectrum 
from which the expected flux in the X-rays and visible can 
be finally obtained. In practice these fluxes are highly vari- 
able (as the gamma-rays) and have been averaged over the 
duration tgo of the burst. 

We did not consider in this paper inverse Compton 
emission models but we now briefly discuss their opti- 
cal properties since they greatly differ from synchrotron 
emission models. When the gamma-rays come from in- 
verse Compton scattering, a fraction a syn ~ 1 — one of 
the total power is emitted at a typical synchrotron energy 
-Ep/r 2 ~ 10 — 100 eV where E p is the peak energy of the 
gamma-ray spectrum. The fraction a syn is fixed by the ra- 
tio Qc/qb through Eq. (17) and (23). Preliminary results 
indicate that for Qc/ob = 1 the emission in the visible of 
a burst of fluence _F 7 and duration tgo — 10 s could be as 
bright as V — 5 — 6 which is already excluded by the limit 
set by the ETC and GROCSE instruments (Krimm, Van- 
derspek & Ricker 1996; Lee et al. 1997). With a c /a B = 100 
the predicted magnitude becomes V = 8 — 9 still within the 
reach of ETC and GROCSE. 

In synchrotron emission models the optical counter- 
part is much weaker. Taking for example the spectrum 



Figure 12. i? p -fluencc relation (thick line) for the burst of Fig. 6a. 
The peak energy is represented in logarithmic scale to show the 
section of linear decline following the moment of maximum count 
rate in the four intensity pulses. The photon fluence corresponds 
to the integrated photon flux in BATSE bands 2 and 3. The thin 
line shows the count rate in the same bands. 

represented in Fig. 8 which corresponds to a burst with 
tgo = 10.55 s we get a V magnitude of 18.4 for F-, ~ 
2 10~ 6 erg.cm~ 2 . A larger fluence would naturally produce 
a brighter counterpart, possibly up to V ~ 15 — 16 for 
Fj <; 2 10~ 5 erg.cm~ 2 . An optical emission in this lumi- 
nosity range should be detectable by the next generation of 
counterpart search instruments, such as the TAROT project 
(Boer 1997). 

In X-rays the calculated fluxes for the BeppoSAX Wide 
Field Camera are in reasonable agreement with the obser- 
vations. In the case of GRB 970228 the X and gamma-ray 
fluences during the first 100 s are F x (2- 10 keV) ~ 1.2 10~ 6 
erg.cm~ 2 and F 7 (40-700 keV) ~ 1.1 10~ 5 erg.cm~ 2 (Costa 
et al. 1997a). For the same gamma-ray fluence the model pre- 
dicts F^ odcl ~ 0.7 - 1.1 10~ 6 erg.cm^ 2 depending on burst 
hardness (the upper and lower limits of F™ odcl respectively 
correspond to bursts with E p — 100 and 350 keV). 

In GRB 970508 the peak X-ray flux and the gamma- 
ray fluence are .Fx (2 - 10 keV) ~ 1.2 10~ 8 erg.cm- 2 .s _1 
and F 7 (50 - 300 keV) ~ 1.1 10~ 6 erg.cm~ 2 (Costa et 
al. 1997b; Kouveliotou et al. 1997). The calculated X-ray 
flux for a burst of comparable hardness, jF™ odcl ~ 8 10~ 9 
erg.cm _2 .s _1 , is close to the observed value while the initial 
X-ray emission predicted by afterglow models appears to 
be an order of magnitude weaker (Sahu et al. 1997b). This 
can be considered as an indication that both the gamma- 
rays and the initial X-ray emission are produced by internal 
shocks. 

The evolution of the afterglow which follows the emis- 
sion from internal shocks can be obtained from a solution 
of the relativistic Sedov problem (Meszaros & Rees 1997b) . 
Such a solution applies when the expanding shell of mass M 
has swept up a mass Mism ~ ™ in the interstellar medium 
which occurs after a deceleration time 

t dec ~ 180 E 5 1 2 /3 n 1 " 1/3 r- 8/3 s , (41) 

in 52 —1 

where E52 is the shell energy in units of -j^- erg.sr , n\ 
the density of the interstellar medium in atom.cm -3 and 
Y2 = jtJq • Most of the energy is radiated at the synchrotron 
frequency of the relativistic electrons 

v s ~ 3 W 16 a 2 e a B Eli 2 t-^ Hz, (42) 

where a c and as are the fractions of the dissipated energy 
which go to the electrons and magnetic field respectively. 
The flux at the synchrotron frequency is given by (Waxman, 
1997a,b). 

F s ~ 6.5 10 _26 y / n7aB-E52-C ) Gpc erg. cm -2 . s -1 . Hz -1 , (43) 
and for v > v B by 

(44) 

where /3 = , q being the exponent of the power-law dis- 
tribution of the accelerated electrons (N(T C ) oc F" 9 ). The 
flux in the visible decreases from a maximum ax ~ !F S fol- 
lowing a power-law of index =|£ the value of (3 (~ 0.6 — 0.8) 
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being obtained from the observations. For GRB 970508 a fit 
of the data also provides Jf™ ~ 6 10~ 2S erg.cm^^s^VHz" 1 
(Sahu et al. 1997b). With the gamma-ray fluence measured 
for this burst and assuming a 5% efficiency for the conver- 
sion of wind kinetic energy to gamma-rays we obtain 

VnTffB ~ 3 10~ 2 . (45) 

If we adopt this value as typical we find that the afterglow 
of the burst (given in example above) with E52 = 2 and 
£>G P c = 2 has a V magnitude of 18.9 before the phase of 
power-law decline. In inverse Compton emission models the 
initial optical counterpart produced by internal shocks is 
considerably brighter than the afterglow but in synchrotron 
emission models the two contributions have comparable 
brightness. However, the optical signal will present a short 
interruption if the duration of the burst is smaller than the 
deceleration time. Conversely, if t w <; idee the two contribu- 
tions overlap and form a continuous signal. In any case, a de- 
tailed photometric follow-up in the optical range beginning 
before the end of the gamma-ray burst would certainly pro- 
vide crucial informations about the emission mechanisms. 



5 CONCLUSIONS 

We have developed a simple model to compute the temporal 
and spectral properties of GRBs under the assumption that 
they originate from internal shocks in a relativistic wind. 
We have not discussed the critical point of how such a wind 
could form but the recent observations of optical counter- 
parts for GRB 970228 and GRB 970508 seem to indicate 
that a relativistic shell was indeed present in these objects. 

The distribution of the Lorentz factor in the wind has 
no reason to be uniform and variations on several timescales 
(down to one millisecond which corresponds to the dynami- 
cal timescale of a relativistic disc orbiting a stellar mass 
black hole) can be expected. Layers of different velocities 
will collide and form internal shocks within the relativistic 
wind and the energy dissipated in these shocks can be emit- 
ted in the form of gamma-rays. We have followed the evo- 
lution of the relativistic wind using an approach where all 
pressure waves are suppressed so that layers only interact 
through direct shocks. We assume that the magnetic field 
reach equipartition in these shocks and that the electron 
Lorentz factor remains close to a constant value r o ~ 10 4 . 
This could be the case if, according to Bykov & Meszaros 
(1996) only a small fraction £ of the electrons is accelerated 
in the shocks and if ( is also approximately proportional to 
the dissipated energy. Gamma-rays are then directly pro- 
duced by synchrotron emission from the relativistic elec- 
trons. 

This procedure allows us to construct synthetic bursts 
whose temporal and spectral properties are compared to the 
observations. We obtain a series of encouraging results: 

1) It is possible to generate a great diversity of profiles 
simply by playing with the initial distribution of the Lorentz 
factor in the wind. 

2) The profile of individual pulses is asymmetric, close 
to a "FRED" shape. 

3) The short time scale variability of the profiles can be 
explained if the Lorentz factor itself varies at the millisecond 
level. 



4) Synthetic spectra can be fitted with Band's formula 
with parameters a, [3 and Eo comparable to those of ob- 
served bursts. 

5) The duration-hardness relation is a natural conse- 
quence of the model since in short bursts dissipation occurs 
closer to the source where the magnetic field and the syn- 
chrotron energy are larger. 

6) Spectral hardness and count rate are correlated du- 
ring burst evolution, the hardness generally preceding the 
count rate. 

7) The pulse width decreases with increasing energy 
following a power law W(E) oc E~ p with p ~ 0.4 for pulses 
of duration 2 - 10 s. 

8) In the declining part of intensity pulses the peak 
energy E p decreases exponentially with photon fluence. 

Some other properties of observed GRBs are however 
not so well reproduced by our model: 

1) The shortest pulses tend to decay faster than they 
rise instead of being symmetric. 

2) Synthetic bursts do not show a global hard to soft 
evolution as frequently as real bursts. 

3) The slope of the E p - fluence relation can differ 
among pulses inside the same burst. 

It is not yet clear whether these difficulties represent 
real problems for our model or are simply a consequence 
of some of the crude assumptions we have made. In parti- 
cular, the evolution of the wind should be followed with a 
detailed hydrodynamical code (Daigne and Mochkovitch, in 
preparation) rather than with the simple method used here. 

A more fundamental issue may be the rather low effi- 
ciency for the conversion of wind kinetic energy to gamma- 
rays. If GRBs result from the coalescence of neutron stars 
the wind cannot be strongly beamed since the merging rate 
is not very much greater than the burst rate. A large fraction 
of the energy released in the coalescence should therefore be 
injected into the wind. 

We believe that the present work has shown that if a 
relativistic wind carrying enough energy can be produced, 
the internal shock model appears as a convincing candidate 
to explain GRBs. Ways to both reach a high efficiency in 
the generation of the wind and to avoid at the same time 
baryonic pollution are among the difficult problems which 
remain to be solved. 
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